V /) 


/as -73 -<?£. 
f'/O/ 


TRITON'S GLOBAL HEAT BUDGET 


f fk 


R . H. Brown , T. V. Johnson, *7. D . Goguen 


Mall Stop 183-501 
Jet Propulsion Laboratory 
Pasadena, CA 91109 


Gerald Schubert 


Earth and Space Sciences Department 
Institute of Geophysics and Planetary Physics 
University of California 
Los Angeles, CA 90024 


Martin N, Ross 


Space Sciences Laboratory 
The Aerospace Corporation 
Los Angeles, CA 90009 


(NASA-CR-138178) TRITON’S GLOBAL HEAT 
BUDGET (JPL) 1& P 


03B 


N9 1-2 30 1 8 


Unci as 

G 3/93 0008401 


PRECEDING PAGE BLANK NOT FILMED 


- 2 - 
Abstract 

Internal heatflow from radioactive decay in Triton's interior along with 
absorbed thermal energy from Neptune total 5 to 20% of the Insolation absorbed by 
Triton, thus comprising a significant fraction of Triton's surface energy balance. 
These additional energy inputs can raise Triton's surface temperature between -0.5 
and 1.5 K above that possible with absorbed sunlight alone, resulting In about a 
factor of -1.5 to 2.5 increase in Triton's basal atmospheric pressure. If Triton's 
internal heatflow is concentrated in some areas, as Is likely, local effects such as 
enhanced sublimation with subsequent modification of albedo could be quite large. 
Furthermore, Indications of recent global albedo change on Triton suggest that 
Triton's surface temperature and pressure may not now be in steady state, further 
suggesting that atmospheric pressure on Triton was as much as 10 times higher in the 
recent past. 
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Voyager's flyby of Neptune revealed that at 38K Triton is hotter than simple 

oh 

thermal models predict (1-3). This conclusion was based A the assumption that 
absorbed insolation is the only significant energy source, and that thermal proper- 
ties common to other icy satellites also apply to Triton (1) . More detailed models 
that include volatile transport also predict surface temperatures on Triton signifi- 
cantly lower than those measured by Voyager, suggesting that Triton has low emis- 
sivity (4,5). Because Triton has a surface composed primarily of solid nitrogen 
(6,7,8), because its atmosphere is composed primarily of nitrogen gas in vapor- 
pressure equilibrium with solid nitrogen on its surface (2,3), and because theoreti- 
cal estimates suggest solid nitrogen has low emissivity (9), this is a currently 
accepted interpretation. 

Despite first-order success in attempts to model volatile transport on Triton, 
there has been some difficulty in accounting for Triton's global albedo distribution 
and the apparent extent of its south polar cap (4,5). Several explanations for the 
discrepancies have been offered, including the supposition that new ice or frost on 
Triton is darker than old ice or frost (4) , but none are very satisfying. Accord- 
ingly, the suggestion advanced here is that difficulties encountered in modeling 
volatile transport on Triton may result in part because significant and previously 
unconsidered, non- solar energy sources are available to Triton's surface. In par- 
ticular, internal heatflow, absorbed radiation from Neptune, and long- time-scale 
changes in bolometric bond albedo may have substantial effects on the global surface 
temperature and pressure of Triton. In addition, anisotropies in Triton's global 
internal heatflow may have substantial effects on local sublimation rates and sur- 
face volatile coverage. 

That internal heatflow on Triton or absorbed thermal radiation from Neptune 
could have measurable effects on Triton's surface temperature and volatile transport 
follows partly because Triton's 38 K temperature is the lowest known of any solar- 
system body (2). This is due to Triton's great distance from the Sun and high 
albedo (1,5), resulting in absorption of only a small fraction of the incident solar 
energy (-1.5 W m ). Also, the ubiquitous solid nitrogen on Triton's surface is 
capable of efficient redistribution of absorbed energy such that Triton should be 
nearly globally isothermal (10). Strong seasonal changes are also now occurring on 
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Triton as it nears a maximum summer solstice in its southern hemisphere (11,12), 
resulting in substantial volatile transport (4,5,10). Further contributing (perhaps 
substantially) to global volatile transport on Triton are several active plumes (1). 

Triton's rather high density of 2.08 gm cm 3 , implies that its mass is -70% 
rock (1), and this in turn implies a large, current source of internal heat from 
radioactive decay. In addition, Triton may have been strongly heated by tidal dis- 
sipation during its hypothetical capture by Neptune (13,14), although it is doubtful 
that tidal heating is presently significant (14) . That Triton has been strongly 
heated is confirmed by Voyager images showing widespread and relatively recent geo- 
logical activity (1) . 

Triton's -15-/ibar surface pressure (2) is sufficient to buffer the temperature 
of the nitrogen- ice portions of Triton's surface such that they are globally isoth- 
ermal (10). This allows a calculation of expected steady-state temperature given 
bolometric emissivity, absorbed insolation and the extent of surface ice or frost 
coverage. A similar calculation was carried out for Triton (1), but the result is 
only approximate because internal heatflow and absorbed thermal radiation from Nep- 
tune were not taken into account. A more complete expression for the global energy 
balance of the Ice or frost on Triton is: 

S 0 (l-A f )a 

+ Ha" + N 0 e f a' - <7e £ T £ a" (1) 

4*-R 2 

where S Q is the total luminosity of the Sun, R is heliocentric distance, a is the 
projected area of frost illuminated by the Sun, a' is the projected area of frost 
illuminated by thermal radiation from Neptune, a'' is the total area of frost, H is 
the global average heatflow, N^ Is the total thermal flux per unit area from Neptune 
at Triton's distance from Neptune, A f is the bolometric bond albedo of the frost, 

is the bolometric emissivity of the frost, T f is the temperature of the 
frost and o is the Stefan- Boltzmann constant. 

More detailed models of frost temperature and volatile transport on Triton 
(4,5) suggest that, in response to seasonal and diurnal variations of insolation, 
some areas on Triton could be stripped of their seasonal coating of frozen nitrogen 
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as summer progresses in Triton's southern hemisphere. At present, however, ground- 
based spectroscopy suggests that the bulk of Triton's illuminated surface is covered 
with nitrogen ice (6,7,8). Thus, for the purposes here, we will assume that Triton 
is completely covered in frozen nitrogen which reduces Eq. 1 to: 


T - 


N 0 * S 0 (l-A) 

— + r 

16*R Z 



( 2 ) 


We drop the subscript f here because we are referring to global properties. Equa- 
tion 2 specifies the global -average temperature on Triton given bolometric bond 
albedo, bolometric emissivity, global -average heatflow and the total thermal lumi- 
nosity of Neptune. The bolometric bond albedo of Triton was derived from Voyager 
images (1,12,15) and from the measurements of the Voyager 2 Photopolarimeter experi- 
ment (16). Although the bolometric bond albedo of Triton derived from the photopo- 
larimeter data disagrees somewhat with that derived from the imaging data, the con- 
sistency of results obtained by several investigators analyzing the imaging data, 
including our own extensive analysis of both groundbased and Voyager imaging data, 
leads us to adopt a value of 0.85 ± 0.05, the approximate average of the published 
imaging results and of our own results (1,12,15). 


Calculation of N 0 , the thermal radiation from Neptune incident on Triton, can 
be made using Neptune's size and effective temperature, and Triton's distance from 
Neptune. The effective blackbody temperature of Neptune is 59.3 ± 2.0 K (3), and, 
ignoring the small change in incident flux due to Triton's finite size, the incident 
thermal flux at any distance from Neptune is given by: 


N 0 - 2aT* 
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for r n « R*. 


(3) 


With r n - 25225 km and ^ - 354600 km (1), N 0 is -3.6 mW m' 2 at Triton's orbit. The 
absorbed fraction of this energy is proportional to Triton's projected area and 
thermal absorptivity (averaged over all wavelengths of interest in Neptune's thermal 
spectrum) which we assume is equal to its thermal emissivity. Finally, for the 
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moment we treat Triton's spherically averaged heatflow as a free parameter, confined 
to the range discussed below. 

Current radiogenic heatflow on Triton depends on the total mass of silicates in 

Triton and the present abundance of long-lived radionuclides. Triton's silicate 

mass fraction probably ranges from 0.65 to 0.75 (1), corresponding to 1.4 to 1.6 x 

10 kg of silicates. For outer planet satellites the most common assumption is 

that silicates have a chondritic composition which produces about 5.45 x 10 W 

kg 1 (17) from radioactive decay. Current estimates of lunar heatflow give a heat- 

“12 “1 

ing rate of 9.20 x 10 W kg (17), which we treat as an upper limit to the 
radiogenic heating in Triton. The resulting range for Triton is then about 3.3 to 
6.6 mW m 2 or a total of 0.75 to 1.5 x 10 11 W. 

The combined contribution of internal heatflow and absorbed thermal radiation 

- 2 

from Neptune is then -4-7 mW m . The relative sizes of different contributions to 
Triton's surface heat budget are listed in Table 1, where the sura of internal heat 
flux and absorbed thermal radiation from Neptune is seen to be 5 to 20% of the inso- 
lation absorbed by Triton. For a bolometric bond albedo of 0.85 ± 0.05, an emis- 
sivity of 0.6 (9) and a heliocentric , distance of 30.2 AU, a temperature increase of 
between 0.5 and 1.5 K is achievable over that resulting from absorbed insolation 
alone. The significance of this is underscored because a 0.5 to 1.5 K increase in 
surface temperature on Triton would increase its atmospheric pressure by a factor of 
1.5 to 2.5! This is shown in Fig. 1, where surface pressure versus emissivity is 
plotted for 3 different values of Triton's bolometric bond albedo; clearly demon- 
strated are the large incremental effects of internal heat flux and absorbed thermal 
energy from Neptune. Thus, Triton seems to be one of the few objects in the solar 
system whose internal heatflow may significantly affect its surface temperature and 
pressure. 

Although the above discussion assumes that Triton's internal heatflow is iso- 
tropic, potentially more interesting is the idea that Triton's heatflow is likely to 
be anisotropic. Support for this comes from images of Triton's equatorial regions 
and low northern latitudes which show definitive evidence of large-scale volcanic 
and tectonic activity (1). If these regions are geologically young, as their lack 
of cratering suggests (1,18), heatflow in these regions could be significantly 
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greater than the global -average value, due possibly to bodies of cryogenic magma 
within a few tens of kilometers of the surface. One way to concentrate heatflow on 
a larger scale would be for convective circulation in Triton's mantle to transport 
internal heat preferentially to the surface in some regions. A permanent, 
kilometer- thick deposit of frozen nitrogen over most of Triton's southern hemisphere 
(as suggested by Voyager images) may act as a heat barrier, causing Triton's inter- 
nal circulation to behave in just such a way. This is conceivable because the ther- 
mal conductivity of nitrogen ice is -50 times less than that of water ice at 38 K, 
thus an extensive cap of frozen nitrogen (such as seems to exist on Triton) could 
act as a barrier to internal heatflow in the regions under the cap material. 
Nevertheless, the plausibility of mantle upwelling in Triton's equatorial regions 
and the role played by a thermally insulating polar cap needs to be assessed by 
numerical models of convection in Triton's mantle. 

If Triton's heatflow is anisotropic, profound effects on local sublimation 
rates and albedos could result, depending on the magnitude of the anisotropy. For 
example, at 60% albedo, areas on Triton near 20° north latitude presently receive 
enough diurnal -mean insolation that a 38 i 1 K temperature can be maintained with 
little or no net sublimation over a Triton day (5). For an emissivity of 0.6, 
between 64 and 79 mW m is radiated to space at 38 ± 1 K. Thus, if local heatflow 
were just a few times the global average, a substantial amount of net sublimation 
would occur at 20° north latitude. If the areas of enhanced heatflow were axisym- 
metric, the condensation front of new frost In the present season would be expected 
to occur substantially north of 20°. In fact, if heatflow in places were -10-15 
times the global average, enough energy would be available to prevent ice condensa- 
tion in any season, regardless of albedo, even for those areas north of the present 
arctic circle on Triton. Thus, areas on Triton showing evidence of enhanced subli- 
mation correlated with past or present geologic activity may be seen in Voyager 
images . 

Up to now we have tacitly assumed that Triton's global bolometric bond albedo 
is constant with time. Differences in Triton's present visual and near- infrared 
reflectance as compared to that in 1976 and 1977 (19,20), however, suggest that 
Triton's global albedo may be variable on a time scale of several years. The data 
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show that Triton was much redder 13-14 years ago than now (1), and thus had a sub- 
stantially lower bolometric bond albedo. There is further evidence for a change in 
Triton's reflectance from infrared observations reported by Cruikshank et al . (21). 


If a change in Triton's bolometric bond albedo has occurred, then it must have 
affected Triton's surface energy budget. Because Triton's atmosphere stores a large 
amount of latent heat, the global surface temperature and pressure on Triton 
require^ a finite time to respond to changes in the global energy budget, thus sug- 
gesting that the temperatures and pressures at Triton's surface recorded by Voyager 
2 may only be transient values. 


To estimate the time scale for the response of Triton's atmosphere to albedo 
changes, we assume global energy balance. Therefore, for a unit area of ice on 
Triton's surface interacting with the atmosphere, energy balance can be written as: 


ueT 


k 

f 



+ S 


(4) 


where h s is the latent heat of sublimation of nitrogen (taken here as 260 J g~ 

[22]), E is the total column mass density of nitrogen gas in the atmosphere, and S 

is the total of absorbed insolation, global heatflow and all other heat sources per 
. 'ft 

unit area averaged over the entire surface of Triton. S can be expressed as: 
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(5) 


Because the global -average sublimation or condensation rate required to maintain a 
temperature of -38 K is -10 gm cm sec , to a very good approximation Triton's 
atmosphere is in hydrostatic equilibrium (23). The surface pressure is then given 
by P — Eg, where g is the mean acceleration of gravity over the atmospheric column 
weighted by the local number density. Equation 4 can then be written as: 


dP _g_ 
dt ” h s 



oe £ 



( 6 ) 


If the atmosphere is always in vapor -pres sure equilibrium with the surface ice, the 
Clausius-Clapeyron relation, P — K exp (-b/T), defines the dependence of basal 
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atmospheric pressure on temperature, and equation 6 becomes (where we again drop the 
subscript f ) : 


dT _ _gf_ 
dt h s Pb 




(7) 


For the mean acceleration of gravity over the atmospheric column on Triton we 
can use as a good approximation the surface value of 78 cm sec calculated from the 
mass of Triton (3) and its radius (1). A detailed calculation of the mean accelera- 
tion of gravity weighted by the density profile of a hydrostatic atmosphere shows 
that the correction for Triton is only -3%. The relation describing the vapor pres- 
sure of nitrogen as a function of temperature is taken from Brown and Ziegler (24). 
The global -average emissivity is assumed constant because temporal variations in 
emissivity are likely to have a much smaller effect on the global energy budget than 
albedo changes* All other parameters are as previously defined except for the time 
history of the global bolometric bond albedo of Triton. We have calculated the 
bolometric bond albedo of Triton for 1976-1977 using the approach of Hillier et ai . 
(15) and the average spectrum of Triton for 1976 and 1977 as published by Cruikshank 
et al . (19) and Bell et ai . (20). From this we find that the bolometric bond albedo 
of Triton in 1976-1977 was -0.74. Because there aren't sufficient data to specify 
the detailed time history of Triton's albedo, we treat 2 simple cases: an instan- 
taneous albedo change and a linear change with time. In Fig. 2 are displayed the 
results of a simulation of the response of Triton's atmospheric pressure to an 
instantaneous albedo change and a linear change over 5 years, where we have included 
a contribution from Triton's internal heatflow and absorbed thermal energy from Nep- 
tune . 


The major aspects of the calculations displayed in Fig. 2 are the initial and 
final atmospheric pressure, and the time scale for relaxation to a new equilibrium 
state after the albedo changes have ceased. Depending on the emissivity of Triton's 
surface, relaxation of the atmospheric pressure in this model can take as long as 15 
years. For lower emissivities , the model suggests that Triton's atmospheric pres- 
sure could have been an order of magnitude higher as recently as 10-15 years ago, 
and that it may be relaxing to a value half or less of the -15 //bars found by 
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Voyager. To effect an increase or decrease of -150 /tbars in Triton's surface pres- 
sure, about 2 gm cm of nitrogen must either evaporate from or condense onto 
Triton's surface globally. This amounts to a change in atmospheric mass and latent 
heat by -4.5 x 10 17 gm and 1.2 x 10 20 J respectively, or an average of -2.5 x 10 11 W 
over a period of 15 years. Thus, it is clear that if albedo changes of the size 
indicated have indeed occurred on Triton, they may still be having a substantial 
effect on Triton's global temperature. 
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Figure Captions 

Fig. 1. Triton's surface pressure as a function of emissivity for three values of 
bolometric bond albedo. The heavy line is for a model that includes solar, internal 
radiogenic, and Neptune energy inputs. The light line is for a model that includes 
only solar input, and the area between the dashed, horizontal lines is the range of 
pressure 15 ± 5 /ibars, encompassing the range of Voyager measurements of Triton's 
surface pressure (2) . Radiogenic heat and heat from Neptune in this model amount to 
6 mW m 2 . 

Fig. 2. Response of Triton's atmospheric pressure to a change in bolometric bond 

albedo from 0.74 to 0.85. The heavy curves are for a linear change in albedo over a 

period of 5 years and the light curves are for an instantaneous change. Curve pairs 

are plotted for several different emissivities , ranging from 0.55 for the top pair 

- 2 

to 0.70 for the bottom pair in steps of 0.05. A combined contribution of 6 mW m 
from internal heatflow and absorbed thermal energy from Neptune has been included, 
chosen so that for an emissivity of 0.55, the final equilibrium pressure is 15 
/ibars . 
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Table 1: Triton Power Comparisons 


Source 

Power 

(lO^W) 

Incident Solar 

85.5 i 

Absorbed Solar 

8.5 

17.0 

Internal Heatflow 

0.75 

- 1.5 

Neptune Thermal 

0.1 

- 0.2 


Note: Incident solar flux was calculated at a distance of 30.2 AU, the absorbed 
solar flux was calculated for a bolometric bond albedo range of 0.8 to 0.9, internal 
heatflow was calculated according to the description in the text and absorbed ther- 
mal energy from Neptune was calculated for the range of emissivity 0.5 to 1.0, 
assuming Triton's thermal absorptivity is equal to its thermal emissivity. 
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